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Abstract 

Theoretical studies of nonequilibrium systems are complicated by the 
lack of a general framework. In this work we first show that a transforma- 
tion introduced by Ao recently (J. Phys. A 37, L25 (2004)) is related to 
previous works of Graham (Z. Physik B 26, 397 (1977)) and Eyink et al. 
(J. Stat. Phys. 83, 385 (1996)), which can also be viewed as the gener- 
alized application of the Helmholtz theorem in vector calculus. We then 
show that systems described by ordinary stochastic differential equations 
with white noise can be mapped to thermostated Hamiltonian systems. A 
steady-state of a dissipative system corresponds to the equilibrium state 
of the corresponding Ifamiltonian system. These results provides a solid 
theoretical ground for corresponding studies on nonequilibrium dynam- 
ics, especially on nonequilibrium steady state. The mapping permits the 
application of established techniques and results for Hamiltonian systems 
to dissipative non-Hamiltonian systems, those for thermodynamic equi- 
librium states to nonequilibrium steady states. We discuss several impli- 
cations of the present work. 



1 Introduction 



Equilibrium statistical mechanics is a fundamental branch of physics with ap- 
plications extending widely through the scientific domain. Based on several 
fundamental laws or hypothesis, the theoretical frames are self-consistent and 
complete. However, most systems in nature are not at equilibrium. Open sys- 
tems, which exhibit free energy dissipation and thus no detailed balance, are 
of special interest. Our understanding of nonequilibrium systems is still primi- 
tive. Due to the lack of detailed balance, the methods used to study nonequi- 
librium processes lack general or systematic approaches and often depend on 
the specific context. Several attempts have been made to develop a unified 
framework and establish a connection between equilibrium and nonequilibrium 
phenomena [TJ [51 [31 IH [3 IHl 13 ■ For example, Oono and Paniconi noticed the 
similarity between an thermodynamic equilibrium state and a nonequilbrium 
steady state, and proposed an operationally feasible phenomenological frame- 
work called "steady-state thermodynamics" |9|. Hatano and Sasa applied the 
formalism to Langevin systems. They and others derived several useful results 

[iniiiiiiis]- 

Continuous efforts have been made but with limited success to map a non- 
Hamiltonian to a Hamiltonian dynamical system [131 [HI [13 [S] • As Morrison 
pointed out, one important reason for the research is that knowing the Hamilto- 
nian structure of a complex system provides guidance for studying the dynamics 

A related topic it is highly debated issue of the existence of potential func- 
tions for dynamic systems without detailed balance [18l[T9] . On the other hand, 
the metaphor of energy landscapes has long be adopted to understand evolu- 
tion of some biological systems [501 [HI [HI [131 [Mj • Recently, Ao has proposed a 
useful way of decomposing a set of stochastic differential equations (SDEs) and 
construct dynamic potentials [25) . Ao and coworkers examined explicitly some 
linear systems for which analytical solutions can be obtained [26]. However, 
the validity of his procedure for nonlinear systems is obscured by the abstract 
presentation. 

For clarity, we we now summarize the content and the main results of this 
work. First, we will demonstrate that for systems possessing steady states, one 
can follow some standard procedure to decompose the effective flux of the corre- 
sponding Fokker-Planck eg nations [2 71 [28] . Ao's procedure turns out to be equiv- 
alent to a special choice of the effective flux. It is essentially a Helmhotz-Hodge 
decomposition for the corresponding stochastic differential equations. The con- 
clusion is not restricted to linear systems. Next we will prove that there exists 
a mapping between stochastic dissipative systems and Hamiltonian systems, 
nonequilibrium steady-state and thermodynamic equilibrium state, respectively. 
To our knowledge this is the first rigorous and general proof of the mapping. 
The mapping reveals common mathematical structures of classes of dynamic 
systems with and without detailed balance. It serves as the theoretical basis 
for 1) the phenomenological steady-state thermodynamics formulations paral- 
lel to equilibrium systems; 2) the simplified derivation of many corresponding 
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relations following standard procedures discussed in any introductory equilib- 
rium statistical physics textbooks; 3) the concept of potential landscapes which 
has been fruitfully applied in biological applications. Unlike that of Oono and 
Paniconi, our discussion is not restricted to steady states and their vicinity, but 
applies to dynamics arbitrarily far from the steady states. Then we will present 
specific examples, and conclude the paper with discussions on the implications 
of this work and potential future studies. 

2 Theory 

Similar to Hatano and Sasa, we will consider processes described by the following 
form of SDEs, 

m 

dxi/dt = Gi(x) + e^gij(x)Cj(t),i = 1, • • • ,ri. (1) 

In general m and n may be different, C,i{t) are temporally uncorrelated, statisti- 
cally independent Gaussian white noise with the averages satisfying < Ci {t)(,j (t) >= 
5ij6{t — r), g(x) is related to the n x n diffusion matrix gg-^ = 2D, where 
the superscript T refers to transpose. For a physical system e is the Boltz- 
mann's constant multiplying temperature. For a non-physical system, one 
can define an effective temperature relating to e. To help the general read- 
ers, let's consider an example reaction series, ^ A B — ^, the chemical 
Langevin equations are fi[/l]/c?t = kg — ki[A] + y/koCi/V^ — \/ki [A]<^2/V^, 
d[B]/dt = ki[A] - k2[B] + y/ki[A]C2/Vfi - Vfe[-B]C3/v^, where il is the sys- 
tem volume [29l [30] . Equations in the form of Eqn. [T] have been applied to 
problems from physics, chemistry, cellular biology, ecology, engineering, finance, 
and many other fields p9 l [3T 1 13§ 1 133] . 

In a recent paper [25] , Ao argued that one can always construct a symmetric 
matrix S and an anti-symmetric one T, and transform Eqn. [Tjinto, 

d-x 

(S + T).- = (S + T).(G(x)-|-6g(x).CW) 

= -Vx0(x)+eg'(x).C(i) (2) 

where is a scalar function corresponding to the potential function in a Hamil- 
tonian system satisfying {d x d(j))ij = {didj — djdijcj) = 0. To uniquely deter- 
mine the transformation matrices S and T, Ao introduced additional constraints 
g'g'"^ = 2S. This is similar to choosing a gauge in electrodynamics (see below). 
Then S and T are determined by 

9 X [(M • G(x)] = 0, (M)-i + (M)-^ = 2D, (3) 

where M = S + T. The above relation also suggests M^-^ = D -f- Q, where Q 
is antisymmetric. One can identify the symmetric and anti-symmetric matrices 
as the friction and transverse matrices, respectively |25|. This transformation 
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has be applied to network studies f331 [33] . On can further introduce the inertial 
term [25] . and transform Eqn. [2] into, 



dx 

It 
dp 



-T(x) . P 
m 

-S(x) . i 



^g'(x)-CW 



(4) 



(5) 



Clearly the above equations reduce to Eqn. [2] in the limit m — so p — )• 0). 
This introduction of artificial inertia is analogous to the widely-used extended 
Hamitonian method used in molecular dynamics simulations [3S1 [371 |3H1 132] • Yin 
and Ao provided detailed analysis of this limit, and obtained the corresponding 
Fokkcr-Planck equation [ID] , 



9tp = V-M-i[eV + V0]p 
which can be further rewritten as, 

dtp = V • [eD • V + e(V • Q'^) + (D + Q) • V0] p 



(6) 



(7) 



To derive the above expression we have used the antisymmetric property of Q, 
and noticed, 



V • (Q • Vp) 



(V • Q) • (Vp) = 

'dQji 



dxi 



V • (V • Q^) 



(8) 



Further studies are needed to examine the general existence and mathemat- 
ical properties of the transformation. In the next sections we will analyze the 
validity and implications of this transformation starting from the Fokker-Planck 
equation corresponding to the SDEs. 



2.1 Flux decomposition 

In this work we will focus on systems for which steady state distributions exist. 
The following discussions follow the procedure used by Graham and by Eyink 
et al closely \27\ [28] , and provide explicit formula alternative to Eqn. [3] for 
constructing Q. Let's consider a generic Fokker-Planck equation, 

dtp{x,t) = V-[-j7p(x,t)+eD-Vp(x,t)] (9) 

Compared to Eqn. 4.1 of Eyink et al. [55], we use a slightly different form by 
proper redefinition of J'. As will be clear in the following development, with this 
redefinition one avoids the requirement that the matrix D has to be divergence 
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free as Eyink et al. assumed. In the following derivations we will not specify the 
form of J', thus the results are independent of the choice of interpretations of the 
SDEs, a subtle issue we will address in the next section. This is another aspect 
in generalizing the procedure of Eyink et al. |28| , where the Ito interpretation is 
adopted. One can define an "entropy" term from the steady state distribution 
S = \n{pss)- Notice that one can decompose the flux term J" in Eqn. IHlinto a 
conservative part TZ and a dissipative part V, 

V = eT) - VS = -L" ■X,n = J -V (10) 

In the above expression we have written in the Onsager form with L'* = eD 
being a symmetric matrix, and X = — V5 being the general force. (In this 
work, we won't distinguish covariant and contravariant tensors in our notations). 
Substituting the expression of pss into Eqn. |9l one has, 

V • [7^exp(5)] = (11) 

Then following Graham, and Eyink et al, one can relate the divergence free 
vector 7?,exp(5) to an antisymmetric matrix J^, 

n.eMS)^Y.^ (12) 

The antisymmetric matrix T can be further expressed in terms of a vector 
potential (see below). So, 

with the new antisymmetric matrix Lf^ = J^y cxp(— 5). Therefore the drift 
term J' can be expressed as, 

J, = ^(Lt^+^^)X,+J2^ (14) 

j ■' 

We essentially reproduced the derivation of Eyink et al but without the need 
for the divergence free approximation on L'* and L". Notice that in the above 
derivations no assumption is made on whether the system is linear. 

Graham discusses how one can construct the matrix J-" and so L" from the 
steady state distribution and current with additional constraints (e.g., satisfying 
homegeneous Maxwell' equation) |2T. The antisymmetric matrix can only 
be determined up to a gauge transformation. This can be seen by noticing that 
j^a' — _|_ T^exp(— 5) also satisfies Eqn. [T^ provided Y^- ^ jXij = 0- Zia and 
Schmittmann derived the explicit form of the anti-symmetric matrix for linear 
Langevin systems [41]. We suggest that the decomposition can be viewed as 
generalization of Helmholtz's theorem, or the fundamental theorem of vector 
calculus, to arbitrary dimensional systems. 



4 



It is well known that the calculus of a SDE is ambiguous depending on the 
choice of interpretation. Consequently, for a given Fokker-Planck equation, one 
can construct different SDEs. While the decomposition of Eqn. [U is formal, 
one can use different definitions of J . The definition of J due to Ito and 
Stratonovich, respectively, 

Jito = G-V-D 

^Stratonovich — (-^^) 

Where G is as defined earlier in Eqn. [1] One may easily define the corresponding 
expressions for other proposed interpretations {e.g. [42l |43l |44] ) . For a special 
choice Jao = G-V-(L°)^, Eqn. [Hgives G = -(L^ + L'^)-X. The subscript of 
i7 indicates that this choice corresponds to the zero-mass limit interpretation of 
Ao [13 . For a system with steady state, Yin and Ao (also see discussions below) 
showed that the steady state distribution is given by Pss(x) oc exp(— 0(x)/e) 
[40] . In the literature this relation is also widely used to define the potential 
[TS] . Substituting this expression into Eqn. [3 and comparing the result with 
Eqn. rni term by term, one can easily identify that S = — L'* = eD, and 
L'' = eQ. 

Since (D + Q) • \7(j) = G and Q has an undetermined choice of gauge, 
a class of SDEs under the zero-mass interpretation can produce a given set 
of steady state distribution and current. A given form of G corresponds to a 
specific choice of the gauge. As it will be apparent in the discussions below, we 
suggest that this freedom of gauge choice is mathematically analogous to the 
gauge freedom in electrodynamics. It is unclear whether the relations given in 
Eqn. [3] are sufficient to determine Q uniquely except for simple cases discussed 
by Kwon et al. [26] . 

2.2 Dependence and independence on Interpretation 

. While Ao has developed the transformation following an unconventional proce- 
dure, the essence is that one may define a new interpretation of the SDEs, which 
has certain appealing mathematical properties as we will discuss below. Given 
the SDEs under one interpretation, there are standard procedures to obtain the 
corresponding equations under another interpretation [45] . Here we suggest a 
procedure to construct the SDEs of the zero-mass limit corresponding to e.g. 
those of Ito interpretation. First construct the Fokker-Planck equation and the 
decomposition [14] of the latter, then the SDEs of the zero-mass interpretation 
are given by the diffusion matrix D and G = — (D + Q)V0. We will show an 
example below. 

In quantum mechanics, one may choose different representations. Analo- 
gously one may interpret a given Fokker-Planck equation with different stochas- 
tic differential equations. This representation ambiguity is related to the gauge 
invariance discussed above. In one aspect, a class of stochastic processes may 
share the same stationary and at least some of the dynamic properties, in an- 
other aspect, for a given stochastic process one may use different forms of mathe- 
matical description depending on how one abstracts the noise source properties. 
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However, one expects that the physicahy measurable quantities should be in- 
dependent of the choice of interpretation, e.5., the values of ensemble averaged 
quantities since the underlying Fokker-Planck equation is the same. This ob- 
servation has practical importance. While the mapping discussed in this work 
allows straightforward theoretical analysis based on knowledges of Hamiltonian 
systems, it is still unclear how one can adopt the zero-mass interpretation in nu- 
merical simulations, and apply the theoretical results. Instead one can perform 
the usual stochastic simulations under the Ito or Stratonovich interpretation 
with the corresponding SDEs. 

2.3 Mapping 

Now let's define a Langrangian, 

£0 = ^mx^ - + X • A(x) 

The Euler-Lagrange equation from variation of Cq, 5Cq — 0, 

d ( dCo\ dCo 



dt \ dxi ) dxi 

results in Eqns. 2] and [5] without the dissipative terms in Eqn. [SJ provided, 

dxj dxi 

Therefore, T is analogous to the electromagnetic tensor [35]. One can also define 
the conjugate momentum, 

d^^Q A , \ 

OXi 

and a corresponding Hamiltonian, 

Ho = P • i - /:o = - ^W)' + 0(x) (16) 
2m 

Next following Zwanzig |47j , we will show that the dissipative terms in Eqn. 
[5] can be replaced by a harmonic bath hamiltonian. First, let's summarize the 
result of Zwanzig. Consider a system described by a set of state variables z = 
(x, p) (for coordinates and corresponding momenta) , and similar bath variables 
Y. Assume that one can define a function 

H(z,Y) =H,(z)+Hf,(z,Y) 

= H.(z) + \[Y- a(z)]^ . K . [Y - a(z)] (17) 



so 



dz/dt = 3 ■V^iH.+Hb) 
dY/dt = J • VY^fc 
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with 

/ 

-/ 

where / is the identity matrix having the dimension of the total coordinates (or 
momenta). Then one can integrate out the bath variables, and obtain (assum- 
ing that the initial bath variables Yp are drawn at random from a canonical 
distribution, i. e., in contact with a heat bath), 



d2t/dt = J • VM^t) + J • W(zJ • F{t) 

+ f drJ- W(zJ • r(r) • W^{zt-r) ■ it-r 
Jo 



where, 



Vza^(z) 

-K •cxp(U -K) • [Yo -a(zo)] 


eT{t-T) 
K ■ exp(iJ • K) 



W(z) 

m 

< F{t) > 
< F(t)F^(r) > 

m 

For our purpose let's choose, 

^b = Y.Y.[ + 2^ljilo'3 - aa(x)/(%/iVa.2^.)) 

Q=l j=l ^ 

Then Eqn. [T8l becomes. 

dxt /dt 



dpt/dt 



V - / rfr7„(r)pt_,/Af + F„(i) 



(18) 



(19) 

(20) 
(21) 



where. 



7a (0 

F„(i) 



TV 



^(Vxaa)^(Vxa) J2 ^ cos(wji) 
aa(xo) \ 



N 

(Vxfla)^ Uaj(O)- 
J = l 
siTl{uJajt) 



cos(a;Qjt) 



+Pa,(0)- 



Let's use a frequency distribution 



0, > 
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where uia is a cutoff frequency. We replace tlie sums over tlie oscillator frequen- 
cies by integrals, 



JV 



and assume that the system momenta do not change significantly within the 
time interval l/wd. Then Eqn. [21] can be approximated as, 



+ v/2^(Vxaa)CaWj (22) 

where 70 — in j (2w^), and Cq. is a random number drawn from a normal Gaussian 
white noise distribution. 

Consider the dissipative terms — + g'(x)C(t) in the special case that 
G = L • X with L a constant matrix, and the diffusion matrix D can be ap- 
proximated as a constant matrix. Then S and g'(x) are constants [55]. This 
corresponds to a linear expansion of a system around a stable stationary point of 
the steady state distribution. In this case the dissipative terms can be replaced 
by a bath Hamiltonian, 



n, = Y 




En I 



(23) 



where n x to is the dimension of g'. The term Va can be geometrically in- 
terpreted as n vectors {vi,--- ,v„} in the to dimensional Euclid space, and 
Sij oc Vj • Vj = |vi||vj| cos(%). 

To prove existence of the bath Hamiltonian in the general case, one needs 
to show that for any given diffusion matrix S = g'g'"^/2, there is at least one 
set of a so S = Va'^ • (Va'^)-^. Notice that the needed number of Oq,, TVq,, can 
be any positive integer. Let's assume that g^^ are at least piecewise analytical, 
and can be expressed as s-th order polynomial around one point Xq, 

ffy- = 5.y (xo) + V5^j(xo) • Ax -I 

Then one can choose 

a" = a?(xo) • Ax + ^Ax^ • a^(xo) • Ax + • • • 

where a" are vectors, and a" are n-th ranked symmetric tensors. Then, 

^"■i ^ ^ E Wikg'ok + {g'jk^g^k +9'ik^9']k) ■ 

^ k 

-f • • • (24) 
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(Va^(VaV).,=^[(a?),(a?) 



(25) 



Let's equate the above two expressions term by term. For the terms in the 
i-th power of Ax, Eqns. [Ml and [25] give n{n + l)/2 x {i + n - l)!/((ri - l)!i!) 
relations, which are hnear equations for the elements of af_)_]^. Exceptions are 
for the zero-th order terms, which give quadratic equations. There are Na x 
(i + n)l/{{n — l)!(i + 1)!) elements of Si°'_^_l need to be determined. Therefore 
the number of variables is no less than the number of constraints provided 
Na > n{n + l)(s + 1)1 {2s + 2n). For the zero-th term, one can always rotate 
the vectors described in Eqn. [23|in the expanded dimensional space, which 
is still an acceptable solution. Especially one can construct a solution so that 
each vector has nonzero projection on each axis. Then the higher order terms 
af with i> \ can be determined (if the number of variables equals the number 
of constraints) or selected (if the number of variables is larger than the number 
of constraints) in sequence. Consequently, in practice one can always relate S 
to a set of scalar functions {a"} to any necessary accuracy. 

Therefore one can map a dissipative system described by Eqns [1] to a ther- 
mostated Hamiltonian system with % = limm-j.o,i.J6->oo(^o + T~ib)- Caution 
should be taken since that the two limits may not commute. This concludes 
our proof. The Hamiltonian mathematically corresponds to Dirac's constrained 
Hamiltonian [48]. It describes a massless particle, coupled to a set of harmonic 
oscillators, moving in a hypothetical scalar potential and vector (magnetic) po- 
tential field expanded in n-dimensions. The magnetic field arises from breaking 
of detailed balance. The steady state distribution is thus given by the Boltz- 
mann distribution of the Hamiltonian system, 

- / ■v^ exp(--H/e) , 

Pss (x, p, Y = . , , ,^ r-TTTT (26) 

J dxdpdY exp(— rt/e) 

or after integrating the momenta, and the bath variables, 

f dpdYexp(-'H/e) , , , , 

= /dxrfprfYexp(-H/e) ^ ^^P^'^/^^ ^^'^ 

where Y represents all the bath variables. Eqn. [^Tjis conjectured by Ao j25j, and 
is consistent with what has been proposed by others O [18]. Here we reach the 
same conclusion naturally, and provide the theoretical basis for further studies. 
It may at first seem surprising that a nonequilibrium steady state with nonzero 
currents can be mapped to an equilibrium state. The zero mass limit is the key. 
This may share interesting analogies with superconductivity. 
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3 Examples 



To further clarify the above abstract discussions, let us consider a simple ana- 
lytically solvable example, 

^^l_^ + V2Ci(t),^=:E-y + V2C2W (28) 

with the diffusion constant matrix set to identity. The matrix Q can be calcu- 
lated from Eqn. [3l 

Consequently, 



-f 










dt \ 


. y J 






-V J 



2 

= (D + Q)"i • G = -(-3a; + y + 2,x~2y + 1)^ 
5 

<t>^\ (2(a; - 1)' + {x- yf + {y- if) + <f>o (30) 

The gauge invariance suggests that the same steady state distribution and cur- 
rent are given by the following stochastic system (under the zero-mass interpre- 
tation), 

AQ.V0 + x/2(^ ^^jjj ^ (31) 
where 

= ( ) "^P^"^) 

with c an arbitrary constant. On the other hand, one cannot make the same 
statement in general under the Ito or Stratonovich interpretation. In their 
analysis based on master equations, Zia and Schmittmann also have a discussion 
highlighting the fact that more than one dynamic system can lead to the same 
steady state distribution and current |41) . 

As another example, let's consider a system adopted by Gomez-Solano et 
al. to test a modified fluctuation-dissipation theorem [49^ . In the experimental 
setup, a bead is confined to a circle, which has a radius as = 1 with proper 
choice of units). The bead experiences a potential sinO and a constant force F. 
The corresponding Langevin equation is 

de/dt = -HQCose + F+V2DCeit) (33) 

with < Ce{t)Ce{t') >— 6{t — t'), and here for simplicity we take e = fc^T = 1. In 
the original work the Ito or Stratonovich interpretation is adopted. The motion 
along the radial direction is highly confined. The model resembles closely to 
rotary molecular motors, where F can be the applied external torque or the 
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noneqiiilibrium driving force from transmembrane electrochemical potential or 
chemical energy (e.g., ATP hydrolysis). Mathematically the model represents 
an even larger class of systems under external force and periodic potential, such 
as linear translational molecular motors like kinesin [SD]. The "steady state" 
distribution (projected to and normalized within a range of 2tt) is 50 , 

plsiS) - r d0'e[(^«(^"(''+^')-^'"'')-^''')/^] (34) 
Jo 

where M is the normalization factor. The quotation mark on "steady state" 
emphasizes that rigorously speaking it is not a steady state in the actual range 
6 = (— cx),cx)). Consequently it is incorrect to directly map the nonequilibrium 
steady state to the equilibrium state of a Hamiltonian system with a potential 
(i?osin0 — F9). Instead let us rewrite Eqn. [33] and include an additional 
equation for the radial coordinate. 



dO/dt = — ^cos( 



dr/dt = —K,{r — 1) 



V2^Ce(i) 



(35) 



with < C,r{t)C,r{t') 6{t ~ t'), < Cr{t)Ce{t') 0, and k » 1. Transform the 
Langevin equations into Cartesian coordinates. 



dx D X xy y 

— = -[K{r - 1) ]- + Ho— - F— 

dt J, J. j,o 



dy 
dt 



D,y X 
r r r 



'2D 
r 

'2D 
r 



{x(r - yCe) 

(yCr + xCe) 



The diffusion matrix and the flux in the Cartesian coordinate is. 



D 



J = 



D 
D 



(36) 



y 



-[«(r-l)-f]f + ffof|-i^ 
-Kr-l)-f]f-Hof^ + F^ 
V D 



(37) 

For this system V • D = 0. One can show that the overall steady state distribu- 
tion is then given by, 



Pss{r, 0) cx rpl^{e) exp 
Then from Eqn. [121 



1 
2D 



K(r- 1)2 



(38) 



Qi2[x,y) = Qi2{xo,yo) 

Pss[r,0) 

Pss{r',e')idy,-dx) 



[J{x\y')-{B-V\npUr',e')] 



(39) 
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where the integration path connects a reference point {xq, yo) and the end point 
(x,y), and Qi2(a;o,2/o) determined with Eqn. [14] at the reference point. Let 
us focus on the case where the motion along the radial direction is confined to 
r = 1{k » 1). Then with the integration path along the unit circle Eqns. 15^ 
reduces to, 

Qi2(cos6',sin6') cx (40) 

Pss \y ) 

One can then construct the corresponding SDEs under the zero-mass inter- 
pretation. Therefore, one can equivalate the original nonequilibrium dynamic 
system to a massless particle moving along a 1-D circular track with poten- 
tial 4> = — lnpg^(0), and a nonuniform magnetic field perpendicular to the 
plane. The magnetic field can be determined from the antisymmetric part of 
(D + Q)^^. This result can be tested experimentally. 



4 Discussions and conclusions 

For systems having a steady state, we demonstrate the relation between the 
transformation introduced by Ao and the decomposition proposed by Graham 
and by Eyink et al. We suggest that the decomposition is an application of the 
generalized Helmholtz theorem in vector field analysis to arbitrary dimensions. 
The theorem states that any vector field in three dimensions satisfying certain 
mild mathematical conditions can be decomposed into a divergence-free portion 
and a rotation-free portion. We further prove that one can map the stochastic 
and generally non-Hamiltonian system into a Hamiltonian system. The mapping 
allows one to study dynamics of a general system with the techniques developed 
for Hamiltonian systems, especially for equilibrium Hamiltonian systems. If 
one can prove the generalized Hemholtz's theorem, which may be related to the 
Hodge decomposition [5T] , then we suggest that the results in this work may hold 
for systems without steady state. The results may even be generalized to discrete 
dynamics [HJ [52]. Further studies are necessary to clarify the mathematical 
structure and to develop numerical algorithms. We hope that this work may 
inspire further dialogue between researchers in nonequilibrium statistical physics 
(especially dynamics of chemical networks) and in other fields such as plasma, 
superconductivity, quantum field theory, and liquid crystal physics. 

In real applications, constructing the exact form of the overall mapping 
Hamiltonian is impractical. Most existing theories on nonequilibrium steady 
states require knowledge of the steady state distributions and currents [571 [2H1 
1321 141] . The procedure of constructing the transformation matrices given in this 
work is of no exception to this requirement. Due to violation of the detailed 
balance in general there is no readily found stationary distribution as for the 
equilibrium systems. In principle Eqn. [3| leads to the Boltzmann-like distri- 
bution. However solving Eqn. [3| may be even more difficult than solving the 
corresponding Fokker-Planck equation directly. The transformation leads to a 
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new zero-mass interpretation of the SDEs. It is unclear how to develop the 
corresponding numerical algorithms. 

The question is then why take the effort of introducing this transformation 
and new interpretation? We want to emphasize that the major result of this 
work is the existence of the mapping itself. Hamihonian systems have special 
mathematical structures and properties. As stated by Morrison [T31 [T7], the 
mapping to a Hamiltonian system allows one to utilize the special properties 
of the Hamiltonian system to guide further understanding of general dynamical 
systems, which may not been so transparent otherwise. Unlike equilibrium sta- 
tistical physics, the lack of a general framework seriously impedes the develop- 
ment of nonequilibrium theories. The current result provides such a framework. 
The mapping allows one to define a partition function, and thus some pseudo- 
thermodynamic quantities for a dissipative system. Then one can readily apply 
well-developed equilibrium results to nonequilibrum systems. The last one or 
two decades have witnessed important progress in nonequilibrium physics stud- 
ies, and several relations concerning the nonequilibrium steady states have been 
derived. With the mapping, we now know that these results, previously de- 
rived through some special techniques and from phenomenological frameworks, 
must exist, simply because there are counterparts for the equilibrium states. 
Let's just mention a few examples. The result of Hatano and Sasa (Eqn. 11 
of [TU]) turns out to be the equilibrium Jarzynski equality [S3]. One can de- 
velop theories on noise-induced phase transitions [Ml [5S] parallel to those of 
thermodynamic phase transitions, noticing that for the former case the noise 
amplitude replaces the temperature (refer to discussions below Eqn. [1]). The 
fluctuation-dissipation relations are originally derived for processes near equi- 
librium [56l [57l [58]. Recent years there are many efforts on generalizing it 
to systems far from equilibrium [T^J [37] [55] [31] . The mapping between a 
nonequilibrium steady state and a thermodynamic equilibrium state makes the 
generalization straightforward [27] [28] [61] . The existence of the mapping also al- 
lows us to derive a Zwanzig-Mori projection formula for general non-Hamiltonian 
dynamics [62' , where defining an invariant measure is a major obstacle in earlier 
attempts ^63i i64j. Biological networks are important nonequilibrium systems. 
They are typically highly inhomogeneous, and full of competing interactions. 
Thus the dynamics of a biological network in some sense resembles that of spin 
glasses but with nonrandom quenching [65] . It is interesting to notice that the 
mapping between a dissipative and a Hamiltonian system suggest possible glassy 
behaviors for relaxation to a nonequilibrium steady-state. 

When a slow classical Hamiltonian system couples through positional co- 
ordinates to a fast and chaotic classical Hamiltonian system, to a first order 
approximation the fast system exerts on the slow system velocity-dependent 
forces described by a symmetric and an antisymmetric tensors |66[ 167] . The 
latter is called " geometric magnetism" , analogous to the " geometric phase" . It 
is not a conincidence that it resembles the theory discussed here. While the 
former emphasizes the dynamic effect, and this (and Ao's) work emphasizes 
the thermodynamic origin, the source of magnetism is due to the fact that the 
slow system is out of equilibrium. For the systems considered here, the fast 
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systems are described by white noises and thus chaotic. Therefore the current 
formahsm, which apphes to non-Hamiltonian systems, can be viewed as a more 
general description of the phenomenon. 

Eqn. [1] contains white noise. For systems subject to colored noise, in prin- 
ciple one can remove them by introducing airxiliary degrees of freedom. The 
mapping introduced by Ao requires existence of noise. However, there is no 
constraint on the noise strength. Therefore one might use the method even 
to extrapolate results of deterministic cases. For systems with spatial depen- 
dence, e.g. the Kardar-Parisi-Zhang equation (S5], in practice one can always 
replace the corresponding stochastic partial differential equations into ordinary 
stochastic differential equations through grid discretization. 

In this work we derived the mapping between classical non-Hamiltonian and 
Hamiltonian systems. The quantum counterpart of Eqn. [1] and its relation to 
the classical Langevin equation has been widely discussed [521 IZHl [ZD HH IZSl 
ITU [751 [751 [77] . It is of theoretical interest to generalize the current work to 
quantum systems especially in the area of quantum dissipation and decoherence 
for a system approaching a nonequilibrium steady state. 

Ao et al. argues that the zero mass limit corresponds to a new interpretation 
of the stochastic differential equations alternative to the usual interpretations 
of Ito and of Stratonovich [751 IZS] • The ambiguity may arise for the following 
reason: even if the diffusion matrix D is independent of x, after the transforma- 
tion, the corresponding matrix T may depend on x nonlinearly. Our mapping to 
a Hamiltonian system suggests that the zero-mass limit interpretation is physi- 
cally self-consistent. As expected, the ambiguity disappears if the Fokker-Planck 
equation is given. 
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